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van der Waals (vdW) interaction plays a fundamental role in the surface-molecules related phe¬ 
nomena. Tuning of the correlated charge fluctuation in the vdW complex is a plausible way to 
modulate the molecules interaction at the atomic surface. We report vdW interaction tunability 
of the graphene-C02 complex by combining the first principle calculations with the vdW exchange 
correlation density functionals and the time evaluation measurements of CO 2 molecules adsorp¬ 
tion/desorption on graphene under an external electric field. The field-dependent charge transfer 
within the complex unveils the controllable tuning of CO 2 from acceptor to donor. Meanwhile the 
configuration of the adsorbed molecule - the equilibrium distance from graphene and O-C-O bonding 
angle - is modified accordingly. The range of electrical tunability is a unique feature for each type 
of molecules. 


van der Waals (vdW) interaction plays an important 
role in the surface-related physics and chemistry [1]. It 
is the decisive factor in the molecular physisorption on a 
surface [2, 3], adhesion between micro-machined surfaces 
in the micro- and nano-mechanical devices [4], friction 
in tribology discipline [5], and the characteristics in the 
vdW heterostructures [6, 7]. At present, surface modifi¬ 
cation, i.e. tailoring or replacing the original materials, is 
a commonly employed method to selectively create spe¬ 
cific vdW complexes for their favorable properties [8, 9]. 

The London dispersion forces involved in a vdW 
complex of molecules adsorbed on a surface is electro¬ 
dynamically correlated to the charge transfer inside the 
complex. This raises the possibility of tuning such vdW 
interaction electrically, e.g. by applying external electri¬ 
cal field, which offers an easy access to the desired fea¬ 
tures while keeping the original materials untouched [10- 
12]. However, we are not aware of any experimental 
reports on vdW interaction showing the electrical tun¬ 
ability. The difficulty of characterizing the vdW complex 
impedes the direct observation. In this report, we demon¬ 
strate the electrical tunability of the vdW interaction for 
graphene-molecule complexes with first principle calcula¬ 
tions combined with the experimental study of the elec¬ 
trically controlled charge transfer for the graphene-C 02 
vdW complex. As a low-noise Dirac fermion material, 
graphene is strongly modulated by the vdW interactions 
with physisorbed molecules upon its transport proper¬ 
ties; the inaccessible interaction is therefore transferred 
into the easily measurable electric signal [13, 14]. The 
transport measurements verify the tunability of vdW in¬ 
teraction between graphene and CO 2 molecules from the 
varied charge transfer from molecules to graphene under 
external fields. 

First, the first principle calculations are utilized to un¬ 
derstand the interaction between graphene and a CO 2 
molecule adsorbed on it [15]. The nonlocal dispersion 


relationship associated with the molecular adsorption on 
graphene is described by using van der Waals density 
functionals reported by Dion [16, 17]. The potential en¬ 
ergy surface scans show the adsorbed CO 2 molecule he 
most stably in the direction perpendicular to the carbon- 
carbon bond of the graphene hexagon, sitting at an equi¬ 
librium distance of 3.2 A with an energy of 170 meV. 
The low adsorption energy reflects a weak interaction 
between the CO 2 molecule and the graphene and its ph¬ 
ysisorption nature. This is also revealed by the density of 
states (DOS) calculations [18-20] as the electronic states 
of CO 2 are superimposed on to these of the graphene 
(Fig. 1A). The weak interaction of the highest occupied 
molecular orbital (HOMO) of CO 2 (the peak showed at 
-9.31159 eV in Fig. 1A) with the graphene states results 
in a small bonding-antibonding splitting of 11.92 meV 
(Fig. IB). The 3-dimensional plots of the wavefuntions 
visually elucidate these abstract features. Figure 1C de¬ 
picts the bonding wavefunction at -9.60327 eV showing a 
weak graphene-CO 2 orbitals overlap while a repulsion is 
noted in its antibonding wavefunction at -9.59135 eV as 
showed in Fig. ID. Furthermore, these interacting CO 2 
states still locate energetically far below the Fermi level 
clarifies the weak charge transfer within the graphene- 
CO 2 complex. Consequently, the van der Waals interac¬ 
tion plays a dominant factor in the graphene-C02 com¬ 
plex interaction 

Now, we unveil the electrical tunability of the 
graphene-C02 molecule vdW complex by carrying out 
first principles calculation with Mulliken population anal¬ 
ysis [22]. The vdW complex is investigated under 
external electric field, E t , applied perpendicularly to 
graphene. In practice, the electric field is usually intro¬ 
duced to graphene surface by applying voltage, Vr to the 
substrate. A finite element simulation is first employed 
to estimate the strength of electric field generated with 
the application of Vt, which shows E t of -0.15 V/A (0.15 
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FIG. 1. Density of States and wavefunctions of a graphene- 
002 complex. (A) DOS plots of graphene-C0 2 complex, 
graphene, and C0 2 molecule. Vertical dash-lines locate the 
their Fermi levels at -4.7840 eV, -4.7788 eV, and -8.6309 eV, 
respectively. The right panel is the zoomed-in details of the 
grey region. (Gaussian broadening energy in DOS plots: 0.2 
eV) (B) Schematic of the bonding-antibonding splitting in 
graphene-C0 2 complex. (C) and (D) visualize bonding (- 
9.60327 eV) and anti-Bonding (-9.59135 eV) Gamma point 
wavefunctions of graphene-C0 2 complex, respectively. Blue 
and red colors denote the positive and negative iso-surface 
(Iso-value: 0.005). 


V/A) is achieved on the graphene surface for Vr of 40 
V (-40 V) (Fig. 2A). The charge transfer from the CO 2 
molecule to graphene, AQ, is calculated with varying E t . 
Figure 2B plots the calculated charge transfer as a func¬ 
tion of electric field. Without applying a tuning field, 
CO 2 molecule acts as a weak acceptor to graphene, re¬ 
ceiving 0.003e charge. The varied A Q is noted as E t is 
introduced, disclosing the nature of tunable charge trans¬ 
fer in the vdW complex. In particular, the CO 2 molecules 
switch their role to donor in the field of 0.15 V/A (Fig. 
2B). The switched polarity of the graphene-C02 molecule 
complex system under the electric fields is also exhibited 
as a change in the electric dipole moments (figs. S10 and 
Sll). 

The predicted tunable charge transfer in the vdW com¬ 
plex, if any, should be measured as varied doping concen¬ 
tration in graphene under the tuning electric fields. In 
order to prove the above theoretical findings, we measure 
the variations of doping concentrations, n^, in graphene 
exposed to CO 2 molecules under various values of Vr- 



FIG. 2. Electrical tuning of van der Waals bonding of C0 2 
molecule to graphene. (A) Finite element simulation show¬ 
ing the direction of electric filed at the surface of 2 A above 
graphene with 40 V applied back gate voltage. (B) First prin¬ 
ciple calculations of charge transfer AQ, bonding distance 
d and bending angle 6 of the adsorbed C0 2 molecule on 
graphene at various tuning electric fields. Inset illustrates 
a absorbed molecule. (C) Charge density distributions in the 
cut planes across the C0 2 at varied electrical fields (tuning 
voltages). From left to right: -0.15 V/A (40 V), 0 V/A (0 V), 
and 0.15 V/A (-40 V). Blue color denotes electron density 
enrichment and red color shows depletion for the iso-values of 
± 0.0001 electrons/A 3 . 

Figure 3 A illustrates the schematic of the device and 
measurement configuration. The gate modulation mea¬ 
surements are carried out along with the time, from 
which the doping concentrations are extracted (Fig. 3B). 
Thermal annealing at ^ 10 -4 Pa is employed in situ to 
regenerate graphene by cleansing the adsorbates before 
each measurement [14]. For Vr ranging from 0 V to 40 V, 
positive increases in the doping concentration nd are ob¬ 
served, which indicate the adsorbed CO 2 molecules act as 
acceptor. Interestingly, for Vt of -40 V, the negative na 
is measured, signifying the CO 2 molecules act as donor. 
These results strongly substantiate the prediction of the 
first principle calculations. 

We can extract the value of the average charge trans¬ 
fer A Q from Fig. 3B with a simple dynamic model 
[24]. Then, doping concentration is expressed as nd(t) = 
N(t)AQ, where N(t) represents the temporal density of 
CO 2 molecules adsorbed on graphene and is described 

N(t) = —(l-e- wt ) + V 0 , (1) 

Pd 

where no is the density of the molecules in the gas phase 
(remains constant after gas injection is stopped), No is 
the initial adsorption density on graphene, p a and pd are 
the adsorption and desorption rates, respectively. By 
fitting Eq. (1) to the measured temporal change in the 
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FIG. 3. Measurements of CO 2 adsorptions at varied tuning 
voltages Vt- (A) Schematic of the measurement configura¬ 
tion. Tuning voltage, Vt, is applied at substrate. (B) Doping 
concentration at various tuning voltages, Vt, from 40 V to - 
40 V. The positive/negative sign of rid denotes hole/electron. 
The solid lines are the fits to Eq.(l). (C) and (D) plot the 
values of averaged change transfer AQ, adsorption rate nop a 
and desorption rate pd, respectively. Solid line in (C) is the 
linear fit. (E) Change of doping concentration after tuning 
voltage switched from -40 V to 40 V. 

adsorption density, the values of A Q are obtained for dif¬ 
ferent Vt- Fig. 3C clearly manifests the charge transfer 
tuned via applied external fields, showing a good agree¬ 
ment with the first principles calculation prediction. Be¬ 
sides, the adsorption and desorption rates are also ex¬ 
tracted from the data fittings (Fig. 3D). The adsorption 
rate enhanced by the increased magnitude Vt, acceler¬ 
ates the adsorption process due to the polarizability of 
the CO 2 molecules. Hence, the similar adsorption rates 
are read for Vt of -40 V and 40 V. In contrast, the des¬ 
orption rate at Vt of -40 V is approximately twice as 
high as compared with that for Vt of 40 V, as a result 
of a weaker vdW interaction between CO 2 molecule and 
graphene at -40 V, reflecting a longer bonding distance 
and lower charge transfer (Fig. 2B). 

It is noted that the application of constant voltage Vt 
is shortly interrupted for an interval of ~ 0.5 min at each 
individual data point in a series of temporal measure¬ 
ment. During the interval, Vt is swept from -40 V to 
40 V in order to obtain gate modulation curves for rid 
extraction. Under these circumstances, the charge trans¬ 
fer from those previously adsorbed molecules is probably 
altered; and molecules freshly adsorbed during the sweep¬ 
ing exhibit varied charge transfer different from that ad¬ 
sorbed under constant Vt - Here we carry out a measure¬ 
ment at an extreme condition to clarify the interruption 
of constant Vt only slightly affecting the analysis of rid at 
constant Vt- Graphene is firstly kept at Vt = —40 V for 


TABLE I. Electrical tunability of the Charge transfer in var¬ 
ious graphene-molecule complexes. Negative and Positive 
signs indicate the role of molecules as acceptor and donor, 
respectively. 


Molecule 

A Q (e) 

- 0.15 v/A 0 v/A 

0.15 V/A 

C0 2 

-0.006 

-0.003 

+0.001 

c 6 h 6 

0.009 

+0.019 

+0.029 

CO 

-0.003 

+0.001 

+0.004 

nh 3 

+0.015 

+0.018 

+0.022 

0 2 

-0.242 

-0.208 

-0.174 


2 hours; a rid of —0.89 x 10 11 cm -2 is measured (Fig. 3E). 
Then Vt is suddenly switched to 40 V; the positive shift 
of rid with a relatively slow rate of —0.4 x 10 10 cm _2 /min 
is observed. Hence, the short interruption of constant Vt 
has only insignificant impact on the extraction of rid for 
Vt- It also suggests that the electrical tuning plays a 
part during the formation of physisorption or the vdW 
complex; that is to say once the complex formed, it is 
stable to some extent. 

In order to clarify the tunable charge transfer, we 
visualize the two-dimensional charge density difference 
distributions in the graphene-CCU complex in Fig. 2C 
[23] under the tuning field of -0.15 V/A (left), zero 
(center), and 0.15 V/A) (right). Here, the charge den¬ 
sity difference of CO 2 adsorption is defined by pAd = 
PC0 2 +Gra ~ PC0 2 ~ PGra, where pC0 2 +Gra , pGra, and 
pco 2 are the charge densities of the graphene with ad¬ 
sorbed CO 2 , the pristine graphene and an isolated CO 2 
molecule, respectively. The tunability of the charge 
transfer can be understood by applying Lorentz force. 
For a zero tuning field, the vdW interaction governs the 
charge transfer between CO 2 and graphene, which can 
be seen from the induced dipole structure between them. 
As the tuning field is switched on, a strong electrostatic 
force induces a comparable order of charge redistribution 
inside the vdW complex; that is to say, the charge trans¬ 
fer is tuned electrically. Under positive (negative) tuning 
fields, more electrons are taken from (given to) the CO 2 
molecule (movie SI). The reversed dipole distribution 
are clearly seen for the negative and positive tuning fields 
(left and right of Fig. 2C). This is associated with the 
change in other parameters for the graphene-CO 2 vdW 
interactions. Particularly, the vdW bonding distance, d, 
is changed from 3.199 A at the natural state to 3.203 
A (3.194 A) for E t of 0.15 V/A (-0.15 V/A). The bond¬ 
ing angle between carbon and oxygen atoms also varies 
accordingly (Fig. 2B). 

We extend the present theoretical investigation to 
other graphene-molecules complexes. In Table.I, we com¬ 
pare the their charge transfer A Q calculated for various 
vdW complexes at the negative, zero, and positive tuning 
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fields [14, 25, 26]. The results indicate that the individ¬ 
ual complexes show different tunabilities. It should also 
be noted that such a unique tunability does not depend 
on the initial conditions of graphene. Thus, additional 
to the conductivity change of graphene stemmed from 
the molecular adsorptions - a commonly employed gas 
sensing mechanism, the observation of the doping tun- 
ability may provide a new approach to the gas molecule 
detection. 

With the complex of graphene and its phsisorbed 
CO 2 molecules, we have discovered the tunability of the 
vdW interaction at external electric fields. The charge 
transfer in a vdW complex is tuned electrically, thereby 
changing other parameters of interaction thereupon. 
Our results entice us to attempt the possibility of 
recognizing the type of molecules with graphene from 
the unique tunable change transfer between them. 
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